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Local time- dependent theory of Einstein - de Haas effect is developed. We begin with microscopic 
interactions and derive dynamical equations that couple elastic deformations with internal twists 
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by the oscillation of the domain wall. Theoretical results are compared with a recent experiment on 
Einstein - de Haas effect in a microcantilever. 
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I. INTRODUCTION 

Einstein - De Haas effect^ consists of the mechanical 
rotation of a freely suspended body, caused by the change 
of its magnetic moment. The latter can be induced by, 
e.g., the applied magnetic field or by rapid warming. The 
Einstein - de Haas effect is a direct consequence of the 
conservation of the total angular momentum (spin + or- 
bital). Consider, e.g., a solid made of N atoms of mag- 
netic moment M = 7jJ, where J = S + L is the operator 
of the total angular momentum of the atom (that includes 
spin S and orbital moment L), 7j = gje/{2mc) is the 
gyromagnetic ratio for J, e < is the charge of electron, 
and gj = l + [2 J( J+l)]-i[J( J+l)+5(5+l)-L(L+l)] is 
the Lande factor. Total angular momentum of the mag- 
net suspended from a string is a sum of N{J) and the 
mechanical orbital moment, due to the rotation of the 
solid. If, for example, the solid, initially non-magnetized 
and at rest, develops a macroscopic magnetic moment 
Al = A^(M) = 7jA/'(J), then the conservation law re- 
quires that N{J) -\- C = 0. This gives C = — AI/7J, that 
is, the solid begins to rotate on being magnetized. 

Experiments on Einstein - de Haas effect and the re- 
lated Barnett effect^ (generation of the magnetic moment 
by mechanical rotation ), performed at the dawn of quan- 
tum physics, provided first measurements of the gyro- 
magnetic ratio for various materials^. Even today the 
Einstein - de Haas method can still provide a more accu- 
rate value of gj as compared to electron spin resonance 
and ferromagnetic resonance methods that require pre- 
cise knowledge of the effective magnetic field inside the 
sample^. Nevertheless fundamental questions about the 
Einstein - de Haas effect remain unanswered. In partic- 
ular, the global conservation of the angular momentum 
does not explain how the angular momentum is actu- 
ally transferred from individual atoms to the whole body. 
This question is clearly related to the magnetic relaxation 
and decoherence at the atomic level. The latter deter- 
mine the width of para- and ferromagnetic resonances, 
as well as functionality of spin-based qubits. Advances 
in manufacturing and measuring of nanomechanical de- 
vices promise to revive interest to the local dynamics of 
Einstein - de Haas effect. 

Our interest to this problem has been ignited by a re- 
cent experiment performed at the NIST laboratory in 



Boulder, Colorado^. In that experiment a 50nm permal- 
loy film was deposited onto a 200/imx20/imx0.6/im can- 
tilever. The cantilever was placed inside a coil that gen- 
erated an ac magnetic field. Oscillation of the cantilever 
was measured by a fiber optic interferometer positioned 
above the tip of the cantilever. When the frequency of 
the ac field matched the resonance frequency of the can- 
tilever the amplitude of the oscillations was about 3nm. 
The data were analyzed within a model that replaced 
the mechanical torque due to change in the magnetiza- 
tion by the effect of the periodic force acting on the fic- 
titious point mass at the free end of the cantilever. Such 
an approximation, while catching some features of the 
phenomenon, is clearly insufficient for the study of the 
microscopic dynamics of the Einstein - de Haas effect. 

In this paper we will develop theoretical framework for 
the description of the dynamics of the Einstein - de Haas 
effect, that we will apply to the problem of magnetic can- 
tilever. To make this problem more transparent we shall 
assume (as is the case for many magnetic solids) that the 
magnetism is of spin origin and can be described either 
by individual spins, S^, localized at the atomic sites 
or by a continuous spin field S(r, t). (Generalization to 
magnetism of spin and/or orbital origin can be obtained 
through a straightforward re-definition of the constants.) 
We shall derive general equations describing the transfer 
of the spin angular momentum to the mechanical angu- 
lar momentum of the body. In the NIST experiment, the 
effect of the ac magnetic field was likely the displacement 
of the domain wall separating two magnetic domains in- 
side the permalloy film. We shall pay special attention 
to this case. The cantilever problem will be solved by 
adding the internal torque due to the motion of the do- 
main wall to the equations of the elastic theory describing 
the motion of the cantilever. The obtained dynamics of 
the cantilever is rather rich and it allows a detailed com- 
parison between theory and experiment. 

General theory of spin-rotation coupling will be stud- 
ied in Section II. Equations of elastic theory with inter- 
nal twists due to the dynamics of spins will be derived in 
Section HI. Mechanical motion of the magnetic cantilever 
will be studied in Section IV. Suggestions for experiments 
will be given in Section V. 
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II. MICROSCOPIC THEORY OF 
SPIN-ROTATION COUPLING 



is known^, this rotation can be equivalently performed 
by the {2S + 1) x {2S + 1) matrix in the spin space, 



Spin-lattice interaction comes from magnetostriction 
and spin-rotation couphng. Only the latter, however, is 
responsible for the Einstein - de Haas effect. The most 
obvious effect of local elastic twists comes from the de- 
pendence of the energy of a spin on its orientation in 
the crystal - magnetic anisotropy. This effect is due to 
spin-orbit interactions and is of relativistic origin. It is 
described by the crystal-field Hamiltonian that can be 
very generally written as 



af3 



af3'yS 
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Here the Greek letters denote Cartesian components of 
a dimensionless spin vector Sj belonging to the site j of 
the crystal lattice. Tensors K^^ ^ j^cx^iS ^ ^^^^ describing 
magnetic anisotropy, are defined in the coordinate frame 



,(1,2,3) 
^3 



' ' ' that is rigidly coupled to the locally-defined crys- 



tal axes, see Fig. [TJ Local rotation of the lattice is per- 




FIG. 1: Rotation of the crystal field due to local elastic twist 
at the location of the spin. 



formed by the (3 x 3) rotation matrix Rj, 



It results in 



(1,2,3) 
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For a small rotation at the site j by an angle 5(j)j^ one 
has 



if 



ra/3 



(4) 



We now notice that due to the rotational invariance 
of Ha^ the rotation of the local frame {e^^\e^^\e^^^) 
is equivalent to the rotation of the vector Sj by the 
same angle in the opposite direction, S R~^S. As 



RjSjRj 



Rj = e 



(5) 



Consequently, 
becomesii^ 



where 



in the presence of rotations, Ra 



RHaR~ 



(6) 



(7) 



In the linear order on Scj)j one obtains 
RHaR-'=Ha+Hr, HR = iY, [Ha, Sj] ■5<j>j . (8) 
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By quantizing 54>j one can apply this Hamiltonian to the 
study of rigid spin clusters and quantum dots^^^^^^^. 

The total spin Hamiltonian Hs may include exchange 
interaction, magnetostriction, Zeeman interaction and 
dipole-dipole interaction. The dipole-dipole interaction is 
usually the weakest one and will not be considered here. 
The magnetostriction is local on spin. Consequently, it 
is transformed by rotations the same way as the crystal 
field. The Zeeman interaction of spins with the external 
magnetic field B, 



(9) 



is not affected by rotations; g being the gyromagnetic 
factor for the spin. Here we take into account that the 
magnetic moment due to spin, M = — ^/i^S, has direc- 
tion opposite to S because of the negative gyromagnetic 
ratio for the electron. Finally, the exchange interaction. 



3 -> 



(10) 



only depends on the local arrangement of spins that is 
not affected by rotations. In the first order on 50^, the 
generalization of Eq. ([8|) is 



E 
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[j?s, S,] • - i J2 [{Hz + He.), Si] • <5</.. 



(11) 

The last two terms appeared in Eq. ([TT]) because Zeeman 
Hamiltonian (|9]) and the exchange Hamiltonian (p!Q|) , that 
are included in Hs^ are independent from local rotations. 
Consequently, one should subtract Hz and ^^ex from Hs 
when computing the effect of rotations. 

Let the total Hamiltonian of the system that incorpo- 
rates all couplings, including interactions with rotations, 
be H. It is clear that the difference between [Hs^^j] 
and [H^ Sj] begins with the terms that are linear on S(pj. 
Thus, in the linear approximation on Scj)j^ we can replace 
i[Hs,Sj] in Eq. with i[H,Sj] = hSj. Working out 
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the commutator with the Zeeman Hamiltonian in Eq. 
(pT]) one obtains 



. (12) 



III. ELASTIC THEORY WITH INTERNAL 
TWISTS DUE TO SPIN-ROTATION COUPLING 

Our approach to Einstein - de Haas effect is based upon 
Eq. ([12]). To apply this equation to long- wave torsional 
deformations of the body we shall write 5cj)j in terms of 
the displacement field of the elastic theory u(r, t)^-, 



(5(/)(r)=ivxu(r), 



(13) 



and replace Sj by the spin density S(r, t). The classical 
energy of the body then becomes 



(14) 



Here He the elastic energy of the body written in 
terms of u(r, t) while Hs = (Hs) includes exchange, 
anisotropy, Zeeman and dipolar energies, magnetostric- 
tion, etc., written in terms of S(r, t) and u(r, t)^^. The 
last term in Eq. (p!l|) follows from equations (p!2|) and 
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fiS + S X (b + bex) ■ (V X u) , (15) 



wher( 



and 
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(16) 



(17) 



The dynamical equation for the displacement field 



(18) 



where (7o,(^ = 5H/5eoc(3 is the stress tensor, ecx(5 = 
dua/dxf^ is the strain tensor, and p is the mass density of 
the material. The stress tensor can be divided into two 
parts, aaf3 = cr^^ 



^(^)--(^) with 



al3 ' 

(M) _ ^{Hs + He) 



and 



SHr 



(19) 



(20) 



Here is the mechanical part of the stress tensor, 

e.g., the part coming from the elastic properties of the 
cantilever plus magnetostriction, while is the part of 
the stress tensor produced by the internal rotations due 
to spins, 

=-^e,;37{^^7+[Sx (b + bex)]^} . (21) 

Notice that, contrary to the symmetric stress tensor 
{o'af3 = o-f3a) uscd by the conventional elastic theory, a^^^ 
is antisymmetric. The immediate consequence of that is 
a non-zero torque. 



dK 



(R) 

af3 



(22) 



acting on the volume element dP'r. Such torques, ne- 
glected by the conventional theory of elasticity, are re- 
sponsible for the Einsten - de Haas effect. 

Equations (p!8|) - ([2T]) allow one to obtain the general 
dynamical equation of the elastic theory that accounts 
for local internal forces due to the dynamics of spins in 
a ferromagnet: 



da 



(M) 

a/3 



J a 



(R) 



where 



fW = _:^v X 
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as + S X (b + bex) 



(23) 



(24) 



Let us check that these equations conserve the total an- 
gular momentum (spin + orbital). Writing the total an- 
gular momentum due to the spins and the crystal as 



J = J d^r [hS + p (r X li)] 



(25) 



one obtains the following equation for the time derivative 
of the a-th component of J, 



J a 



= J d^r ^hSa + eaf3-frf3 {pil^) 
= J d^r ^hSa + eaf3-frf3\/s(T, 



(M) 
7(5 



hSf3 + (S X b);3 - ef3^sIer^S^\/e'^r^Ss\ } ,(26) 



where we have used equations (j23j) , ([24j) and (p!6|) . If one 
prohibits transfer of spin angular momentum through the 
surface, integration by parts with account of the symme- 
try of 0-^^^ and 7^/3 gives 



with 



j = K(^)+K(«), 



dA^ 



(M) 



(27) 



(28) 
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and 



(29) 



Here K^^^ is the external mechanical torque applied to 
the surface of the body A, while K^^^ is the volume spin 
torque due to the external magnetic field. Thus, in ac- 
cordance with our expectation, when external forces are 
absent, Eq. ([23|) conserves the total angular momentum, 
j = 0. 

If the spin-lattice interaction was absent there would 
be no deformation induced by the dynamics of spins. 
This condition provides another check of the validity of 
equations ([23|) and ([24|) . In the absence of dissipation the 
spin field satisfies^^ 



-S X h 



'eff 



(30) 



where beff = 6H/6S is the effective field that can be pre- 
sented as beff = b + bex + b^ Here b' is determined 
by the the spin-lattice coupling. Its main part is usu- 
ally the anisotropy field, b^ = SHa/^^- Equations ([23|) 
and (|3Q|) then show that Zeeman and exchange interac- 
tions alone do not provide any force on the body. This 
is in accordance with the fact that spins should couple 
to the lattice in order to produce such a force. Dissi- 
pation can be incorporated into the problem by adding 
standard damping terms to the elastic equation ([23|) and 
Landau-Lifshitz equation (|3Q|). 



IV. DYNAMICS OF MAGNETIC CANTILEVER 

Among many problems involving internal forces due to 
spins, Eq. ([23|) can be used for computation of the elastic 
motion of a magnetic cantilever. For example, in the case 
of the motion of a domain wall inside the cantilever, one 
substitutes the known domain wall solution for S(r, t) 
into the right-hand side of Eq. ([23|) . while the left-hand 
side follows from the elastic theory of cantilever in the 
absence of spins^^. Sudden increase of the external mag- 
netic field should result in the domain wall sweeping the 
cantilever, thus providing a source of deformation during 
a finite time. Application of a harmonic ac magnetic field, 
as in the NIST experiment, should lead to the oscillation 
of the position of the domain wall inside the cantilever. 

The geometry of the problem is shown in Fig. O The 
cantilever of length L, parallel to the F-direction, is mag- 
netized in the X-direction. The ^ = end of the can- 
tilever is attached to the holder while the y = L end 
is free. We are interested in small displacements of the 
cantilever in the Z-direction, Uz{y^ t), caused by the time- 
dependent external magnetic field. Vectors B and S are 
assumed to lie in the XY plane. The latter property of 
the magnetization is common for thin films. It is easy 
to see that in this case the terms proportional to S x b 
and S X hex in the right-hand side of Eq. ([23|) give zero 
contribution to the Z-component of the elastic equation. 
Adding the term proportional to S to the right-hand side 



Domain Wall 



FIG. 2: Geometry considered in the paper. 



of the conventional equation of motion for a cantilever—, 
one obtains from Eq. ([23|) 



P- 



h^E d^u. 



h d d 



dt^ 12(1 -cr2) dy^ 2dydt 



sAy.t). (31) 



where h is the thickness of the cantilever in the Z- 
direction, E is the Young's modulus, and a is the Poisson 
coefficient, — l<cr<l/2. 

If the magnetization of the cantilever was rotating uni- 
formly in space, then according to Eq. (|3T]) the force from 
spins would only act on the free end of the cantilever 
where the magnetization has a discontinuity. However, 
for a soft magnetic material like permalloy, deposited on 
a cantilever that is large compared to the dimensions of 
a monodomain particle, the change in the magnetization 
should occur through the motion of a domain wall. For 
that reason we shall describe the magnetic state of the 
cantilever by two domains separated by the domain wall 
dit y = yo{t). Then Sx{y^t) is given by the domain wall 
solution centered at = yo{t)^ 

S.{y,t) = Sd^[y-yo(t)] . (32) 

In the absence of the dc magnetic field this gives 



df^ 12(1-0-2) %4 



(33) 



For the dissipative motion of the domain wall the speed 
of the wall is proportional to the field. When the ac mag- 
netic field, B = Bo cos{ut)^ is applied in the A-direction, 
one has yo{t) = ^o(O) + asm{u;t)^ where a < L is the 
amplitude of the oscillations around yo{0). The domain 
wall is given by 5'dw(^, t) = SqF [y — yo{t)] where So is a 
constant spin density and F changes from —1 to +1 as 
one crosses the wall. Note the connection of Sq to the 
magnetization, Mq = qiibSq. 

It is convenient to switch to dimensionless variables, 



-, t = tv, 



~ Y 12p(l - (j2)L4 ' 

(34) 

where L is the length of the cantilever and u determines 
the scale of the eigenfrequencies of its vertical oscillations. 



Uz{y,t). In terms of these variables Eq. becomes 



(9 9 f dyo\ O F 

= —e 



dt J df ' 



where 



hSo hSo /12p(1-(t2) 



(35) 



(36) 



is a dimensionless parameter that does not depend on the 
length of the cantilever L. By order of magnitude, e ~ 
h/{Msh)^ where M ^ p/So is the mass of the material 
per spin 1/2 and s ~ y/E/p is the speed of sound. It is 
easy to see that e is a small parameter that can hardly 
exceed 0.01 even for the smallest cantilevers. 

For the given function yo{t) that describes the motion 
of the domain wall, Eq. ([35|) has to be solved with the 
following boundary conditions: 

du 

Uz = 0, = at ^ = , 
dy 



dy^ 



at y = l. (37) 



The first two conditions correspond to the absence of 
displacement and the absence of bending of the cantilever 
at the fixed end, while the last two conditions correspond 
to the absence of torque and force, respectively, at the 
free end^^. 

For the free oscillations of the cantilever, e = 0, one 
writes 

Uz{y, t) = u{y) cos{ujt) . (38) 
Substitution into Eq. ([35|) with e = then gives 



dy^ 

The general solution of this equation is 



(39) 



u{y) = A cos{hzy) + B sm{hzy) + C cosh.{tzy) + D smh.{hzy) , 

(40) 

where A^B^C^D are constants of integration. With the 
help of the first, second, and fourth boundary conditions 
(|37|) one obtains 

^ , ^ ^ ^ sin n — sinh n , , , 
C = -A, D = -B, B = —A. (41) 

cos K: + cosh K: 

Substitution of these expressions into Eq. (|4Q|) gives up 
to a normalization factor 

u{y) = {cosK + cosh/^) [cos{hzy) — cosh.{hzy)] 
-\-{smK — sinhA^) [sin(/^^) — sinh(Ai:^)] . (42) 

The third of the boundary conditions (|37|) provides equa- 
tion. 



cos hi cosh -\- 1 = , 



(43) 



n=l 
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FIG. 3: Profiles of the oscillating cantilever at different mo- 
ments of time for n = 1, 2, 3. 



for the frequencies of the normal modes of the cantilever, 
On = i^n (measured in the units of u of Eq. (|34j) ). Fun- 
damental (minimal) frequency is Coi ^ 3.516. The next 
two frequencies are 002 ~ 22.03 and Os ^ 61.70. The pro- 
files of the oscillations of the cantilever for three normal 
modes (n = 1, 2, 3) are shown in Fig. [3l 

We shall now turn to the forced oscillations of the can- 
tilever due to motion of the domain wall. We first neglect 
dissipation and write for the displacement 



(44) 



where Rm{t) are functions of time to be determined and 
Um{y) are normalized eigenfunctions (|42|) of the free can- 
tilever that correspond to eigenvalues tZm given by Eq. 
(031), 



/" 

^0 



dyum{y)u„{y) = 5„ 



(45) 



Substitution of Eq. (011) into Eq. gives 
'd'^Rr 



^\ dp 



u^lRm ) u^{y) = -e ( S > (46) 



dt J dy'^ 



where we have used Eq. ([39]) . Multiplying both parts of 
this equation by Un{y) and integrating over y from to 1 
with account of Eq. (05]) . one obtains linear second-order 
differential equation for Rn{t), 



d Rfi 
~dP 



dyo 



oiRn = -M ^ ) / dy—un{y) . (47) 



When the width of the domain wall is small compared to 
the length of the cantilever, the first derivative of F can 
be replaced by the ^-function. 



dF 
dy 



25[y-U^- 



(48) 



6 



In this case, integrating by parts in the right-hand side 
of Eq. (|47|) , one obtains 



d Rfi 



2e 



di 



dUn ' 



-2e^^in [yo{t)] . 



(49) 



Dissipation can be included into the problem by adding 
the first time derivative of Rn to this equation. This 
results in a conventional problem of damped oscillations 
induced by a periodic force: 



dRn 



d^ Rn 



dt^ Qn dt 



ulRn = -2ej-Un[y^{t)] . (50) 



Here Qn is the quality factor of the oscillations of the 
cantilever at the eigenfrequency udn- 

The most interesting case is when the position of the 
domain wall, 



yoit) = 6 + asin(a;t) , 



(51) 



oscillates at a frequency uo that is close to one of the 
resonant frequencies of the cantilever o;^ = To solve 
Eq. ([5Q|) we write Un and Rn as Fourier series, 



Rn{t)= J2 



in) ikOt 



Substitution into Eq. ([5Q|) gives 



ikoJOJri 



(52) 
(53) 




0.0 0.1 0.2 0.3 



0.4 a 0.5 



FIG. 4: Color on line: The dependence of P^^\a,b) on the 
amplitude of the oscillations of the domain wall whose equi- 
librium position is in the middle of the cantilever: a) k — 1, 
h) k = 2. 



where 



1 

2^70 



27r 



u)^'^ (a, b) = ^ I d^e-'^^Un (6 + a sin ^) . (54) 



Writing u^^-' as u)^^ = | exp [^7^"'^] one obtains the 
following expressions for the amplitude |r["^^| and phase 



Xn) _ I ,(n)| 



of the k-th harmonic of Rn{t): 



(»)i 



2efc(I>|4"^| 



(55) 



where 



2w„(l) 



(58) 



The dependence of P^^^ on d and b for various k and 
n is illustrated in Fig. [4] and Fig. [5l Note the non- 
monotonic dependence of the amplitude of the free end 
on the equilibrium position of the domain wall for n = 2 
and n = 3. It is due to the profile of the normal modes 
of the cantilever shown in Fig. [3l When equilibrium 
position of the wall coincides with the antinode, the effect 
of the oscillation of the wall on the cantilever is minimal. 



An) 



in) 

Tk 



arctan 



kcJCJn 



(56) 



According to Eq. ([55|) . resonances occur at frequencies 
uj = uon/k that are independent of damping. Ai uo = 
ujn/k the maximum displacement of the free end of the 
cantilever is given by 



(57) 



V. SUGGESTIONS FOR EXPERIMENT 

Expressions derived in this paper provide the frame- 
work for theoretical analysis of the experimental data on 
Einstein - de Haas effect. To illustrate applications of 
the theory we have derived rigorous formulas for the me- 
chanical motion of a magnetic cantilever, induced by the 
motion of a domain wall when the cantilever is placed in 
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FIG. 5: Color on line: The dependence of P^^\a,b) on the 
equilibrium position of the domain wall for n = 1,2,3 and 
/c = 1,2. 



the ac magnetic field. In our theory we assumed that the 
entire volume of the cantilever was magnetic. The for- 
mulas can be easily adjusted, however, to the situation 
when the magnetic layer has thickness ph < as was the 
case in the NIST experiment. In this case the strength of 
the source in the right hand side of Eq. (|35|) reduces by 
the factor p. Consequently, one should replace e in the 



above formulas with Cp = pe < e. 

Accurate comparison between theory and experiment 
requires precise knowledge of the mechanism by which 
the magnetic moment is changing. If it is due to the mo- 
tion of the domain wall, as we believe was the case in the 
NIST experiment!^, then one needs to know the initial 
equilibrium position b and the amplitude of the oscilla- 
tions of the wall a. The parameter b can be controlled 
by a weak dc magnetic field, while a can be controlled by 
the amplitude of the ac field. It is also desirable to excite 
various harmonics cOn/k and to identify the fundamental 
frequency uoi . This would allow one to obtain the value of 
the parameter u in equations and ([36|) . If the mag- 
netization and the ^-factor are known, the ratio So/p in 
the first of Eq. (|36|) can be computed with good accuracy. 
The precision with which the parameter e can be deter- 
mined will then depend on the knowledge of the length of 
the cantilever L. Alternatively, e can be extracted from 
experiment and used to obtain the spin density So- If 
the magnetization Mq is known this would allow one to 
obtain the gyromagnetic factor g = Mq/ {iibSq). 

In the NIST experiment the fraction of the magnetic 
material p was close to 1/12 while the dimensions of the 
cantilever were L = 2 X IQ-'^m, = 6 X 10"'^m. This 
gives ep ~ 10~^. If a is comparable to L, then according 
to Eq. ([57|) , the deflection of the free end of the can- 
tilever at the fundamental frequency uo = uoi must be of 
order CpQiL. The observed deflection in the nanometer 
range then corresponds to Qi ~ 100. We should notice 
in this connection that the effect could be stronger for a 
cantilever with a higher quality factor. As a matter of 
fact the quality factors as high as 10, 000 have been re- 
ported for microcantilevers^^. For a cantilever of length 
L = 0.2mm such a high quality factor would allow the 
deflection of the free end due to Einstein - de Haas effect 
as high as a few tens of a micrometer. 
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